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ABSTRACT 

From  the  equations  of  motion  for  two  charged  fluids  with 
scalar  pressures  it  is  shown  that  there  exists  a  steady  shock 
profile  if  there  is  friction  between  the  two  fluids  proportional 
to  the  difference  between  their  velocities. 
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Introduction 

Shocks  as  strict  mathematical  jump  discontinuities  occur  in 
magnetohydrodynamics  as  well  as  in  classical  fluid  mechanics  if  the 
motion  is  governed  by  conservation  laws  without  friction  termsc   If 
friction  terms  are  present,  or  introduced  in  the  differential  equations, 
the  shock  discontinuities  are  resolved  into  zones  of  rapid  transition, 
and  the  behavior  of  the  flow  in  these  transition  zones,  the  "shock 
profile"  can  be  studied,  much  as  ordinary  boundary  layer  flow  is  studied, 
by  greatly  magnifying  the  width  of  the  shock  layer  and  considering  the 
flow  in  the  shock  layer  as  a  continuous  solution  of  the  corresponding 
differential  equations  which  joins  the  states  on  the  two  sides  of  the 
layer.   These  states  in  turn  are  connected  with  each  other  by  the 
classical  shock  relations,  which  are  expressions  of  the  conservation 
laws.   Such  an  analysis  of  the  shock  profile  is  necessary  to  justify 
using  discontinuous  shocks  as  physically  realistic  models 0   In  magneto- 
hydrodynamics  the  problem  is  complicated  by  the  fact  that  in  many 
physically  interesting  cases  the  geometrical  dimensions  within  which 
the  flow  takes  place  are  small  compared  with  the  mean  free  path,, 

The  present  paper  describes  such  a  profile  analysis  under  a  trial 
assumption  of  a  friction  term  which  may  be  too  much  of  an  ad  hoc 
idealization  but  allows  rather  complete  discussion  and  analysis,,   If 
there  is  no  friction  then  there  is  no  shock  but  there  are  non-trivial 
steady  state  solutions.   The  theory  is  carried  out  for  a  two  phase 
continuum  fluid  (positively-charged  ions  and  negatively  charged 
electrons)  with  scalar  pressures  in  a  one-dimensional  steady  situa- 
tion, in  which  all  quantities  depend  only  on  the  space  variable  x 
and  a  shock  zone  connects  two  given  states  at  -  co  and  +ooa   For  results 
using  a  modified  one  phase  fluid  see  [1],  [2]  and  their  bibliographies. 
For  other  results  using  a  two-phase  fluid  see  [J]» 

The  method  consists  of  reducing  the  problem  to  the  solution  of  one 
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ordinary  nonlinear  differential  equation  of  first  order  with  singu- 
larities. 

The  author  gratefully  acknowledges  the  advice  and  assistance 
of  R,  Courant.   The  computations  were  carried  out  on  the  Univac  by 
Ho  Goertzel  and  the  preliminary  hand  computations  were  made  by 
Carole  Siegel  and  Milton  Vassell. 
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lo   Equations  of  Motion 

We  shall  consider  the  equations  of  motion  expressing  the  usual 
conservation  laws  for  the  two  fluids  together  with  Maxwell's  equa- 
tions 0   We  assume  that  the  fluid  is  neutral  which  is  an  extremely 
accurate  approximation  to  Maxwell's  divergence  equation  for  the 
electric  field.   This  yields  a  system  that  is  simple  to  handle  and 
physically  reasonable.   We  also  assume  that  there  is  an  interaction 
between  the  two  fluids  which  may  be  interpreted  as  a  friction  or 
resistance.   Since  the  flow  depends  on  x  only,  these  equations  (see 
Grad  [3]  for  their  general  form)  in  this  special  case  reduce  to: 


(1) 


d!(V^   =  ° 


(Conservation  of  Mass) 


(2)  £.   (p   u^-u1   +   p11) 

v     '  dx     Kr  r  r  r    ' 


s=l,2 
(Conservation  of   Momentum) 


(3)        dx    (erur"   +   uJPr1}   =   E*(^rur)    +     &r        (Conservation  of  Energy) 


(k) 


(6) 


d!*1    -    o 


#•  E2  =  A  E3 
dx  dx 


=      0 


(5)  (curl  H)1   =   J1   =  J-  q   u1 

T    r  r    f 


j 


(Maxwell's    equations) 


where 


(7) 


qr     =      Wr 


Here,  i  =  1,2,3  denotes  the  component  in  the  x,y,  or  z  direction! 
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r  =  1,2  denotes  the  ion  and  electron  fluid  respectively,,   pp  is 
density,  u1  is  velocity  and  P^  is  the  pressure  which  we  take  to  be 
scalar, 

(8)  p£J  =  5ijpr  ' 

J1  is  current  density,  E   (or  E)  is  electric  field,  H   (or  H)  is 

magnetic  field,  u.  is  the  permeability  constant,  q^  is  charge  density 

and  y  is  the  constant  ratio  of  charge  density  to  mass  density 
r 

(  6  /m  ),  a   is  the  resistance  coefficient  with  a^2  =  a21  =  a  >  0„ 
For  the  present  we  shall  permit  a  to  be  a  function  of  any  of  the 
dependent  variables  except  the  pressures0   If  it  should  depend  on 
the  pressures  then  it  may  only  depend  on  p^  +  p2»   This  approach  can 

be  modified  for  a  depending  separately  on  p.  or  p?°      er  is  total 

—         12  12 

energy  per  unit  volume   ep  =  ep  +  fPrur  =  cvpr  +  ^  ppur  . 

The  energy  sources  £/   satisfy  the  equation  21  tsv  ~   °>  that  is, 

the  total  energy  equation  is  unaffected  by  the  internal  resistance 

forces  a   (u  -  u  ).   To  determine  the  electric  field  or  the  distri- 
rs  r    s 

bution  of  temperature  one  would  need  more  knowledge  of  o_o   This 

will  be  done  for  the  case  £,     =0, 

r 

We  consider  only  the  case  where 

(9)  H1  =  0, 

i0e0,  there  is  no  component  of  magnetic  field  normal  to  the  shock. 
Finally,  from  the  assumption  that  the  total  charge  density  is 
zero 

(10)  qi  +  q2  =  Jjl   YrPp  =  0„ 


The  boundary  conditions  at  -co  are  that  the  velocities  and  fields 
should  be  constant  and  we  may  choose  coordinates  so  that 
u2,  =  u3  =  0,  H3  =  0,  E2  =  0. 


2 

A 

A  H 


->-l 


u 
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2,   Integration 

Several  of  the  equations  can  be  integrated  singly  or  in  combina- 
tions „   From  (1),  (5)  and  (10)  we  find 

■~Y  Y 

(11)    pu1  =  M,   p  =  Pl  +  p2,   plBS_-p,   p2=:-_p, 


with. 
(12) 


111 

U   =  U,  =  Up 


the  speed  of  each  fluid  in  the  x-direction;   M  is  constant.   By 
adding  the  two  momentum  equations  (2)  in  the  y-  and  z-directions 
respectively  and  using  (12),  (8)  and  (10),  we  find 


_d_ 

dx 


i\   1 

Pxur)u 


q  (u  X  uH)1  = 


*r  r 


(  H  vV  *^H 


(i=2,3) 


By  (9),  (10)  and  (12)  the  right  hand  side  vanishes,  and  hence, 
from  the  conditions  at  -oo, 


(13)        Z  Pr>%   =   © 


r  r 


(i  =  2,3). 


Taking  appropriate  linear  combinations  of  (2),  we  have  by 
(8),  (12),  (13), 


_d_  /  Ti  1 
dx 


(j  u  )  =  -r1r2PE    +  y-lY2pm-u  e1<5 


1     1      (Y-L-Yp)2   , 
+  YnYoPuu  e.  ,HJ  +  a  — = — - —  jx 


PY-^        ' 


where 


e^H^  =  H3 


=  -H' 


for  i  =  2, 
for  i  =  3. 
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That    is,    by   (11),    (5)    and    (1+), 


2 

dH" 


YjYg 


(11,)         ^  =  ^f  ^(uA  £  ff>)    +   u.u  .  -^-2  -g-       ^ 


(15)    E3 . .  _uL  d (ui  d  h2)  _    ^2  +  a(Yri2)    <4^  fs! 

YtYoM  dx    dx       ^  2  2      M2    dx 

' 1 ' 2  YXY2      M 


whe 


re,  by  (6),  FT   and  E  are  constants,.   Without  loss  of  generality 


2  3 

we  may  take  E  =  0.   Then,  from  (1^.),  we  see  that  H   =  0  is  a  solution 

and  satisfies  the  constant  conditions  at  -co.   Thus, 

(16)  J2  =   E2  s  H3   =   0. 

We  have  reduced  the  differential  equations  of  motion  (1-10)  to 

1      2 
(15)  where  u  and  H   are  the  dependent  variables,  and  in  addition 

there  are  two  momentum  equations  in  the  x-direction  (2)  and  two 

energy  equations  (3).   We  shall  use  only  the  total  momentum  and 

total  energy  equations  to  find  the  shock  profile,,   The  two  remaining 

equations  could  be  used  to  determine  the  charge  separation  field  E 

and  the  individual  pressures. 

The  total  momentum  equation,  found  by  adding  the  two  equations 

of  (2),  integrates  to 

(17)  Mu1  +  p  +  I  u-(H2)2   =   P  , 

where  p  =  p,  +  p~,  and  P  is  the  constant  of  integration.   The  total 
energy  equation  (3)  integrates  to 

(18)  (e  +  pju1  =  E3H2  +   I 
where  I  is  the  constant  of  integration. 
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Here 

e  =  cv(Pl+p2)  +jl  Prur*ur 

=  cvp  +  I(p1+p2)(u1)2  +  |  Z   (Pl(^)2  +  P2(^)2) 

1      2  P 
=  cvp  +  -  Mu  -  gj^r-  (dF"} 

by  (12),  (11),  (13)  and  (111).   We  may  therefore  write  ( 18 )  as 

(19)    [(oT+l)(P  -  |^(H2)2  -  Mu1)  +  ^Mu1  -  gjJi—Cg-)  Ju1  -  EV  =  I 

The  problem  has  been  reduced  to  solving  (15)  and  (19)  where 

-i  o  3 

u  ,H  are  the  dependent  variables  and  M,E  , I,P  are  constants  of 

integration.*"'  The  total  pressure  p  is  found  from  (17),  PpJ  up,  J 
from  (5),  (ID,  (12),  (13). 


;  2 

*Note  that  a  =  a(uj,  pp,  H1,  px+p2)  =  a(M,P,u1,H2,  ^-  ) 
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3.   Conditions  on  Shock  Solution 

We  shall  look  for  a  solution  of  (1.5)  and  (19)  subject  to  the 
condition  at  x  =  - od  that  the  state  is  constant  and  show  that  the 
state  to  which  one  is  led  as  x  ->  +oo  is  the  corresponding  state  behind 
a  shock  provided  the  friction  coefficient  a  >  0. 

The  state  at  x  =  -oo  satisfies 


P  =  PQ,   P  =  PQ,   H2  =  HQ  and  u1  =  uQ  >  j/^H2  +  (l+c;1)po/po 

dH2 
(shock  speed  is  greater  than  characteristic  speed).   Also  -r- —  =  0  for 

2 
x  =  -oo  since  H  ->  H  .   From  the  standard  theory  of  magneto-hydro- 
dynamics (see,  for  example,  Friedrichs  and  Kranzer  [$])   we  know  that 
these  constants  determine  the  flow  behind  a  shock.   An  appropriate 
profile  theory  must  lead  us  to  this  flow  at  x  =  +oo. 

The  constants  of  integration  are  found  from  (12),  ( 1$ ) ,  (17), 
(19): 

M=pu,   P  =  Mu  +p  +  iuH  =  p  „u  (1  +  §.  ) 

Ko  o'         o   *o    2r  o   *o  o      p.2 

1    o      -d       i  r>  2+p  (c  +1) 

(20)   I  =  (cT+l)pouo  ♦  lMU2  -  E\   =  lp0«3(i  ♦    °2  Y     ) 

3  3-1-2 

EJ  =  -  [luK      =   -  p  u;?H  A 
r   o  o      o  o  o 

/ p" 2~~  2 

where  A  is  the  Alfven  Mach  number,  V  p^u  /uH   ,  and  (3   =  2p  /(xH   . 

/   o  o   o         o     oo 

We  introduce  some  convenient  dimensionless  variables: 

1^2/ .  a(yn-Yo)2  u        1/2 

(        Vo/-YlY2^o  (-T^M)^"         » 

(21) 

„2  1 

TT        _        H  U 

H  -  —  ■       u  =  —  • 

o  o 
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C=/-YlY2^po  X 


Prom  (15)  and  (21),  we  obtain  a  pair  of  differential  equations  for  Y) 
and  H, 


(22) 


u|5  =-l  +  uH-  kurj 


dH 


an     yi 
udT  =  ? 


Here  u  is  a  function  of  H  and  Y)    given  through  the  quadratic  equation 
(19)  which  reduces  by  (20),  since  c   =  — — r-  ,  to 

(23)    A2(y+Du2  -  yu[2A2+1+P0-H2]  +  [A2-  Y|  2-2  (H-l )  ]  ( Y-D  +yPq  =  0. 
The  appropriate  root  is  determined  by  the  initial  conditions: 

(214.)        u-1  =  H-l  =7=0,      £,     =   -cd. 

Since  u  is  continuous  we  may  switch  to  the  other  root  only 
where  the  two  roots  are  equal,  that  is,  by  (17),  (20)  and  (23),  where 

(25)  u2   =  ^   " 

PUo 

The  state  behind  the  shock  can  be  easily  shown  from  the  magneto- 
hydrodynamic  shock  conditions  to  satisfy 

(26)  uH  =   1,         ^  =  0. 

We  note  that  (22)  defines  a  direction  field  on  a  two-sheeted 
surface  with  coordinates  yi  and  Ho   On  each  surface  u  is  defined  by 
the  appropriate  root  of  (23)  and  the  two  sheets  are  joined  where  (25) 
holdSo   The  direction  field  has  singularities  where  (26)  holds.   Solu- 
tions for  which  u  changes  sign  are  physically  meaningless,,   Likewise  u 
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is  bounded  since  by  (17),  the  pressure  p  would  become  negative  for 
large  u„   The  condition  p  >  0  reduces  by  (17),  (20)  and  (21)  to 

(27)       2A2(u-l)   <   1-H2   +   pQ  . 

Finally,  in  general,  no  physically  meaningful  solution  crosses 
the  curve  C  where  (25)  holds.   For,  on  C,  the  two  roots  of  (2I4.)  are 
equal  and  hence  at  any  point  on  C,  ■*■£>   and  -r?   are  the  same  on  both  sheets, 
It  can  be  shown  that,  with  a  few  exceptional  cases,  solution  curves 
issuing  from  >?  =  0,  H  =  1  do  not  meet  C  tangent ially  and  therefore  to 
go  from  the  first  sheet  to  the  second  dH  or  dri     must  change  sign. 
Hence  d£  changes  sign  and  the  solution  is  not  single-valued  in  Z,   (or  x). 
Therefore 

(27a)         u2   >  YP/pu2  . 

We  therefore  will  never  find  a  solution  in  which  the  velocity 
is   less   than  gas  sound  speed.   In  particular  the  only  possible  shock 
profiles  found  by  this  method  will  be  those  for  which  the  state  behind 
the  shock  satisfies  this  condition.   This  limits  the  range  of  shock 
strength  in  terms  of  Alfven  speed  that  we  can  consider.   But  infinitely 
strong  changes  in  pressure  are  not  excluded  since  the  pressure  ahead, 
or  6   can  be  as  small  as  we  like  and  the  pressure  behind  is  finite. 
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i\.B      Nature  of  Singularities 

The  singularities  of  (22)  occur  where  (26)  holds.   By  (23)  there 
are  three  such  points.   One  point  corresponds  to  the  state  ahead  of  the 
shock,  one  to  a  state  behind  the  shock  and  one  to  a  negative  shock 
which  we  do  not  need  to  consider.   The  behavior  near  the  singular 
points  is  given  by 

(28)      ^  =  Re  ae^  +  ...,     H  =  1+Re  be^  +  ... 
where  £  satisfies 


r~i 

(29)    K   =  "  |±  J    \   +  u"2(uH)H 

evaluated  at  the  singular  point. 

For  the  singular  point  corresponding  to  the  state  ahead  of  the 
shock,  we  find  from  (2l|), 


K   .  .  \   +  jM  +   2A2  -  2  -  x£ 


k   is  the  value  of  k  ahead  of  the  shock, 
o 


Since  the  shock  speed  relative  to  the  gas  is  greater  than  the 

2 
characteristic  speed,  we  find  2A   -  2  -  yB  >  0.   Hence  £  is  real,  and 

this  point  is  a  saddle  points 

For  H  =  H..  say,  n  =  0,    corresponding  to  the  state  behind  the 

shock,  we  find  using  (20),  (21),  (23)  and  (2t|), 


„  .    kl  +  /  kl  +  2A1  -  2  -  Igl 

*--W  T+    (2A|.YPi) 

where  k,  and  u,  are  the  values  of  k  and  u  behind  the  shock,  A,  is 


1U  - 


the  Alfven  Mach  number  behind  the  shock, 


2     p(u  ) 
1     ^FF 


=   A  u 


2„3 


and  p,   = 


*5— >?  =  -2-  6   behind  the  shock. 

/„^^    p   'o 


Since  the  speed  of  the  shock  is  less  than  the  characteristic 
speed  behind  it 


2k\      <   2   +   YPr 

P        P 
By  (27a)  only  shocks  with  U-,  >  yp/pu   may  be  considered 


2AX   >   yP1. 

Hence  the  singular  point  corresponding  to  the  state  behind  a  shock 
is  a  focal  point  since  k,  >  0  and  small  as  we  shall  see  in  Section  9. 

n  Direction  Field  in  r\ ,  H-l  Plane 


H-l 


kn  t   0 


H-l 
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5o   Existence  depending  on  k^  f   0 

It  can  be  shown  that  if  u  >  yp/pu   there  will  always  be  a  con- 
tinuous solution  from  one  singular  point  to  the  other  with  d£  increas- 
ing if  k:  /  0. 

In  the  important  limiting  case  k,  =  0  the  focal  point  becomes  a 
center  and  there  is  no  solution  that  enters  this  point,   We  see  from 
(22)  and  (23)  that  the  solution  in  the  (Y[,H) -plane  is  symmetric" 
about   ^  =  0.   The  only  solution  issuing  from  the  initial  state  is 
one  which  leads  back  to  it  (solitary  wave)  see  Figures  la,b.   Thus 
there  does  not  exist  a  shock  solution.   We  note,  however,  the  existence 
of  steady  state  solutions  which  are  periodic  and,  of  course,  do  not 
satisfy  the  initial  conditions. 

When  the  shock  solution  exists  it  will  look  in  the  ( \  ,H) -plane 
as  in  Figure  2a  and,  in  the  £,H-plane,  the  shock  profile  is  as  indi- 
cated in  Figure  2b. 

Finally,  there  is  at  most  one  curve  leaving  the  saddle  point  (0) 
and  entering  the  focal  point  (1).   For,  if  two  curves  leave  the  initial 
point  (0)  with  £  increasing  and  enter  the  end  point  (1),  the  state 
behind  the  shock,  then  so  does  a  curve  leaving  (0)  with  £,   decreasing. 
Therefore  close  to  (1)  these  three  curves  circle  (1)  and  enter  the 
focal  point.   They  intersect  >7  =  0  with  >?  increasing  Infinitely  often 
and  they  never  intersect  each  other.   Therefore  d  r)  /d£  changes  sign 
infinitely  often  on  ^  =  0  or  by  (22)   -1  +  uH   vanishes  infinitely 
often  on   ^  =  0  which  contradicts  (2i|). 


Similar  results  (under  slightly  different  assumptions)  have  been  obtained 
by  A.  Banos,  by  Davis,  Lust,  and  Schluter,  by  Ho  Grad,  and  by  M. 
Mittelman. 
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6.   Wave  Length  and  Shock  Width 

The  wave  length  of  the  oscillations  as  £,   ->oo  is  given  by 

The  wave  length  in  the  x-space  by  (21)  is 

-1/2 

Note  that  the  characteristic  length  X  =  (-Y-iYpP  P-)      Is  independent 

of  the  initial  state  except  for  the  density. 

We  find  next  the  order  of  magnitude  of  the  damping  coefficient 

using   either  Kolodner   [6]  or  Spitzer   [I4.].    By  [I4.]  equations 

2-1).,  p.  18,  with  (2)  we  see  that 

P.e  =  a(u*  -  u|)   =   a(v*  -  v*) 

Comparing  [2]  equation  (2-13)  with  (2-9)  of  this  paper  with  Z  =  1, 
n.  =  n  and  setting  Spitzer' s  n  =  N 

1        e  or- 

N  =  (_2_)2   ?ie   =  (_c_)2  a 
ene    lvi"vel     eni 

Here  c  is  speed  of  light,  e  is  charge  in  e.s.u. ,  n  =  no.  density, 
and  N  is  given  in  e.m.u.  by  [1+]  page  86; 


N  =  1.29  x  10^  i^n. 


Then 


en      2                  in.  /ra.m 

k«  (--£)    -^-  (-£),  x  -iJLJL 

c  p  u       m  / 

ro   o        e  e/p   p, 


or 
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where    e    is   given   in  e.s.u.,    |j,   in  e.m.u.;    or,    since  m     «   m.  , 

y  I        ^  v    „    e2     NX 

/n.|j.  ceo 

These  quantities  are  both  independent  of  the  prevailing  magnetic 

field.   Taking  n,  =  lO^/c.c.,  rn  =  9.035  x  102   gm.   e  =  1^77  x  10"10 

x  e 

abs.  e.s.u.  \i   =  1.26  x  10   henries/meter,  T  =  10   electron  volts, 
we  find  X  =  .017  cm.  and  k  =  1.15  x  10  u"   sec/cm,  with  u  in  cm/sec. 

Let  bX  by  the  distance  from  P,  to  Pp  where  ?,  is  the  point  with 
the  lowest  value  of  x  for  which  u-1  =  ,10(u,-l)  and  Pp  be  a  point  such 
that  I u-u, I  <  .10u  for  x  >  x(Pp).   Then  by  numerical  computation  it 
is  indicated,  for  the  case  with  zero  pressure  ahead  that  with  k  =  l/lOO, 
b  will  be  about  200.   The  details  of  the  first  cycles  are  in  Figure  2b. 

If  b  is  about  200  the  shock  width  is  about  I4.  cm.  which  is  a 

reasonable  width  for  many  Sherwood  applications.   Based  roughly  on  the 

<         / 
above  examples,  this  corresponds  to  u  ro    ICr  cm/sec,  that  is,  to  a 

slow  shock.   In  other  words,  only  if  the  magnetic  field  and  pressure 

ahead  of  the  shock  are  small  will  the  shock  profile  found  here  be 

significant.   On  the  other  hand  as  the  speed  of  the  shock  becomes 

large  the  shock  width  bX  rapidly  becomes  too  large  for  Sherwood 

applications.   This  failure  will  occur  in  any  theory  base  on  collisions. 

It  is  necessary  to  study  the  existence  of  such  fast  shocks  using  the 

Boltzmann  equations  without  collisions  instead  of  a  continuous  theory. 

These  restrictions  do  not  prevent  us  from  getting  large  pressure 

ratios  in  the  fluid  by  having  sufficiently  low  fluid  pressure  compared 

to  magnetic  pressure  (low  B)  ahead  of  the  shock.   In  this  sense  strong 

shocks  are  admitted.   The  shocks  have  the  character  that  their  speed 

is  not  too  different  from  the  characteristic  speed  of  the  medium  ahead, 

and  in  this  sense  the  shocks  considered  are  weak. 
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